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de consumo de nutrientes y crecimiento de un bloom de 
microalgas en un fotobiorreactor High Rate Algae Pond 
(HRAP). Tecnología y Ciencias del Agua, 6(3), 15-31. 


En el presente trabajo se cultivó un bloom de microalgas en 
un fotobiorreactor piloto tipo high rate algae pond (HRAP), 
con un volumen de 450 1. Éste se operó en discontinuo, con 
un medio sintético y bajo condiciones ambientales de luz 
y temperatura. El objetivo era utilizar un modelo cinético 
que predijera las velocidades de crecimiento de biomasa, 
consumo de nitrógeno y fósforo, y la evolución del contenido 
de estos dos nutrientes en la biomasa. Los resultados 
parecen indicar que el contenido en nitrógeno y fósforo 
se redujo en más de un 80% en el agua residual sintética y 
que efectivamente el modelo resultó útil para predecir el 
consumo de nutrientes por las microalgas, lo que junto con 
los datos experimentales del contenido en nutrientes de la 
biomasa indicó la existencia de otros procesos que compiten 
con las microalgas en la eliminación de los nutrientes del 
fotobiorreactor. 


Palabras clave: High rate algae pond, microalgas, nitrógeno, 
fósforo. 


Introducción 


Uno de los problemas actuales de las masas de 
aguas naturales es la eutrofización. La eutrofiza- 
ción se define como “el aumento de nutrientes 
en el agua, especialmente de los compuestos de 
nitrógeno o de fósforo, que provoca un crecimien- 
to acelerado de algas y especies vegetales supe- 
riores, con el resultado de trastornos no deseados 
en el equilibrio entre organismos presentes en el 
agua y en la calidad del agua a la que afecta” 
(Directiva 91/271/CE). Los principales nutrien- 
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tes que conducen a la eutrofización de las aguas 
naturales son las formas de nitrógeno y fósforo 
(N orgánico, P orgánico, NH de NO, NO, y 
PO?) presentes en los efluentes secundarios de 
las plantas de tratamiento de las aguas residuales 
(Ruiz-Marin, Mendoza-Espinosa, €: Stephenson, 
2010). 

Con el fin de minimizar este problema, 
la legislación europea establece una serie de 
requisitos para los vertidos procedentes de las 
instalaciones de tratamiento de aguas residuales 
urbanas realizados en zonas sensibles cuyas 
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aguas sean eutróficas o tengan tendencia a serlo 
en un futuro próximo (Directiva 91/271/CE). 
Los valores límites más restrictivos establecidos 
para fósforo y nitrógeno total son de 1 y 10 mg 
1, respectivamente (Directiva 98/15/CE). 

Las estaciones convencionales de tratamien- 
to de aguas residuales están diseñadas para eli- 
minar de manera eficiente la materia orgánica y 
los sólidos en suspensión, pero no para eliminar 
el contenido en nitrógeno y fósforo total (Abdel- 
Raouf, Al-Homaidan, « Ibraeheem, 2012). Por lo 
tanto, para cumplir con los niveles establecidos 
por la legislación vigente en las zonas sensibles, 
se están incorporando nuevas tecnologías a los 
diagramas de flujos de las estaciones depurado- 
ras de aguas residuales (EDAR) (Arbib, 2012). 
El cultivo de microalgas se propone como un 
tratamiento terciario eficaz para eliminar el 
nitrógeno y el fósforo presentes en las aguas 
residuales (Oswald € Gotaas, 1957). Los meca- 
nismos empleados por las microalgas para la 
eliminación de nutrientes pueden deberse a pro- 
ceso bióticos, relacionados con la asimilación de 
nutrientes, o a procesos abióticos asociados con 
stripping del amonio y precipitación del fósforo 
a valores de pH elevados. Entre las ventajas del 
uso de microalgas destacan el abaratamiento de 
los costos; la eliminación simultánea del nitró- 
geno y fósforo sin necesidad de adicionar una 
fuente externa de carbono; la obtención de un 
efluente oxigenado, beneficiando así el medio 
receptor, y la posibilidad de su implantación 
aguas abajo de la estación depuradora sin ne- 
cesidad de remodelar la planta de tratamiento. 
Además, esta vía también puede ser explotada 
para la producción de biomasa con un alto 
valor como materia prima en la fabricación de 
biocombustibles de tercera generación, al igual 
que para la biofijación de dióxido de carbono y 
su retirada de la atmósfera (Arbib, 2012). 

Los sistemas de cultivo de algas en suspen- 
sión se clasifican por lo general en sistemas 
abiertos y cerrados (Tredici, 2004). Los sistemas 
de cultivo abierto, que se han venido empleando 
desde los años de 1950, se pueden clasificar en 
aguas superficiales naturales, como estanques, 
lagunas y lagos, y estanques artificiales. Dentro 


de estos últimos, los denominados “raceway 
ponds” o “High Rate Algae Ponds (HRAP)” 
son los más empleados. Por otra parte, los sis- 
temas cerrados son reactores transparentes, de 
plástico o vidrio, y se clasifican en placa plana, 
fotobiorreactores tubulares y sistemas de bolsas 
(Cai, Park, €z Li, 2012). Los fotobiorreactores 
HRAP son estanques abiertos, poco profundos 
(15-30 cm de nivel de agua), en los que las al- 
gas y los nutrientes son impulsados por unas 
paletas rotatorias que giran alrededor de un eje 
horizontal. Estos sistemas se caracterizan por 
ser baratos, fáciles de construir y operar (Arbib, 
2012; Cai et al., 2012). 

Con respecto a la eliminación del nitrógeno 
y fósforo del agua residual, lo importante no es 
la velocidad a la que crecen las microalgas sino 
la velocidad a las que las microalgas consumen 
el nitrógeno y fósforo, ya que no todos los 
nutrientes consumidos por las microalgas son 
empleados para crecer. Las microalgas tienden 
a consumir nitrógeno y fósforo en exceso, pro- 
ceso conocido como “luxury uptake” (Powell, 
Shilton, Chisti, € Pratt, 2008). 

El “Phototreatment model (PhBTY es un 
modelo cinético que combina los procesos de 
crecimiento de la biomasa y asimilación de nu- 
trientes y como consecuencia permite predecir 
la evolución del contenido de nutrientes en la 
biomasa en experimentos en discontinuo. Los 
experimentos llevados a cabo para desarrollar el 
PhBT se realizaron a pequeña escala (utilizando 
frascos Pyrex de borosilicato de dos litros como 
fotobiorreactores), y en condiciones controladas 
de luz y temperatura en el laboratorio. Debido 
al pequeño tamaño de los fotobiorreactores, las 
microalgas no se cosecharon, por lo que no se 
comprobó de modo experimental si la cantidad 
de nitrógeno y fósforo presente en las microal- 
gas coincidía con las cantidades estimadas por 
el modelo PhBT (Ruiz et al., 2012a). 

En el presente trabajo se pretende reproducir 
el procedimiento cinético establecido por el 
PhBT, usando un fotobiorreactor piloto tipo 
HRAP, operando en discontinuo bajo condi- 
ciones ambientales de luz y temperatura para 
predecir las velocidades de crecimiento, el con- 
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sumo de nutrientes, y el contenido de nitrógeno 
y fósforo presente en la biomasa. 


Material y métodos 
Microorganismos y medio de cultivo 


Para este experimento se utilizó un bloom de mi- 
croalgas, cuya especie dominante era Coelastrum 
sp., obtenido de los fotobiorreactores HRAP 
presentes en la planta piloto de la EDAR “El 
Torno”, situada en Chiclana de la Frontera, Cá- 
diz (36* 25' 38.15” N, 6” 9 23.9” O), España. Para 
el cultivo de las microalgas se empleó medio 
COMBO (Kilham, Kreeger, Lynn, Goulden, éz 
Herrera, 1998). 


Fotobiorreactor 


Para el desarrollo del experimento se empleó 
un fotobiorreactor piloto tipo HRAP. Se utilizó 
un HRAP de fibra de vidrio, con un volumen 
de 450 l, operando con una lámina de agua de 
aproximadamente 30 cm, y una superficie total 
de 1.93 m? (2 525 mm de longitud, 750 mm de 
anchura y 450 mm de profundidad) (figura 
1). Lo que supone una relación superficie- 
volumen (S /V) de 3.62 m*. El cultivo se mezcla 
mecánicamente con cuatro paletas que giran en 
torno a un eje horizontal movidos por un motor 
a una velocidad óptima entre 5 y 7 rpm (Dodd, 
1986). 

La rueda de paletas está asentada sobre una 
depresión en el fondo del fotobiorreactor que 
sirve para reducir el flujo de retorno. Además, 
las paredes del HRAP están curvadas excén- 
tricamente en sus extremos más alejados de la 
rueda de paletas para acelerar el flujo (Arbib et 
al., 2012). 


Diseño experimental 


El HRAP se puso a funcionar en las instalaciones 
habilitadas en el exterior del Centro Andaluz de 
Ciencia y Tecnología Marinas (Cacytmar). 

El fotobiorreactor se inoculó con el bloom 
de microalgas hasta obtener una concentración 
inicial de entre 0.05 y 0.07 g SS 1". 
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Figura 1. Planos y diagramas de flujos del fotobiorreactor 
HRAP, representados en mm. 


El experimento se mantuvo en discontinuo 
hasta alcanzarse el estado estacionario; esto es, 
cuando la concentración de biomasa algal se 
mantuvo constante durante tres días consecu- 
tivos. 


Procedimiento analítico 


Se tomaron muestras todos los días y de manera 
aproximada a la misma hora. Los análisis rea- 
lizados fueron concentración y composición de 
biomasa, concentración de nutrientes, y como 
parámetros de control se determinaron pH, 
intensidad lumínica y temperatura. 

La biomasa en peso seco se determinó en 
forma de sólidos en suspensión mediante el mé- 
todo normalizado de los sólidos en suspensión 
secados en estufa de 103 a 105 *C (APHA-2540 
D), para ello se emplearon filtros de fibra de 
vidrio (GF/F CAT núm. 1825-047, Whatman'"). 
Ésta también se determinó de forma indirecta 
por densidad óptica a 680 nm. 

Cada vez que la concentración de biomasa 
en el fotobiorreactor se veía incrementada en 
más de 0.1 g SS 1”, se centrifugaba un volumen 
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de cultivo adecuado para obtener la cantidad 
de biomasa necesaria para realizar, previo 
secado por liofilización, su análisis elemental 
de carbono, nitrógeno, hidrógeno, azufre (ana- 
lizador CHNS Leco * 932) y fósforo (digestor 
microondas modelo Speed Wave de Berghof y 
espectrómetro de emisión atómica de plasma 
IRIS Intrepid de Thermo Elemental). Cada uno 
de los elementos analizados se determinó por 
triplicado; en el caso del fósforo, por duplicado. 

La concentración de nitrógeno y fósforo total 
se determinó mediante el método de oxidación 
propuesto por Kóthe y Bitsch (1992). Para la 
oxidación previa de la muestra se añadió una 
microcuchara de reactivo OXISOLV? (Merck 
KGaA, Darmstadt, Alemania) a 10 ml de 
muestra; se incubó a 105 *C durante 60 minu- 
tos y se dejó enfriar a temperatura ambiente. 
Posteriormente se prosiguió con el análisis del 
nitrógeno y del fósforo en forma de nitratos (N- 
NO), y fosfatos (P-PO,*) mediante los métodos 
colorimétricos de Múller y Weidemann (1995), 
y APHA 4500-P-E, respectivamente. 

La concentración de nitrógeno en forma de 
nitratos (N-NO,,), amonio (N-NH,?) y nitritos 
(N-NO,,) se midieron al inicio y final del expe- 
rimento, conforme a los métodos colorimétricos 
de Miller y Weidemann (1995), APHA 4500- 
NH, D y APHA 4500-NO, B, respectivamente. 
La concentración de fósforo en forma de fosfatos 
(P-PO,?*) se midió de modo colorimétrico de 
acuerdo con el método normalizado APHA 
4500-P-E. 

El valor de pH se midió de forma electroquí- 
mica (pH-metro GLP 32, CRISON9). De manera 
puntual se tomaron datos de pH cada minuto 
durante 24 horas mediante una sonda de pH 
conectada a un data-logger de adquisición de 
datos, para poder estudiar cómo varía el pH en 
relación con el ciclo día-noche en función de la 
fotosíntesis. 

También se tomaron datos de la radiación 
fotosintéticamente activa (PAR) y de la tempe- 
ratura tanto externa como interna del fotobio- 
rreactor cada cinco minutos; la adquisición de 
los datos se llevó a cabo mediante un data-logger. 


Análisis estadístico 


El análisis de los datos experimentales de 
evolución temporal de las concentraciones de 
biomasa, nitrógeno y fósforo se realizó confor- 
me al modelo cinético PABT (Ruiz ef al., 2012a), 
basado a su vez en los modelos de crecimiento 
de biomasa microalgal y consumo de nutrientes 
en aguas residuales, propuesto por Verhulst 
(Verhulst, 1838) y Romero (Quiroga, Perales, 
Romero, € Sales, 1999), respectivamente. Para 
hacer el modelo se emplearon técnicas de re- 
gresión no lineal mediante la minimización del 
cuadrado de los residuos empleando la herra- 
mienta Solver de Microsoft Excel 2007. 


Nomenclaturas 

X: concentración de microalgas en un ins- 
tante f. 

X. concentración inicial de microalgas. 

X : concentración máxima de microalgas que 


el sistema alcanzaría en batch. 

u(X): máxima tasa de crecimiento específico. 

E tiempo. 

Pp: productividad volumétrica del fotobio- 
rreactor Operando en discontinuo. 

Py. productividad volumétrica del fotobio- 
rreactor Operando en continuo. 

Pa: productividad superficial del fotobio- 
rreactor Operando en continuo. 

O: tiempo de retención hidráulico (TRH) al 
que opera el fotobiorreactor. 

Mrsp: Profundidad de la lámina de agua del 


fotobiorreactor. 

S: concentración de substrato soluble en un 
instante f. 

S: concentración de substrato inicial (nu- 


trientes en el medio de cultivo). 
concentración de substrato no asimilable. 
X_/Y: concentración de nutrientes en el inóculo. 
St: concentración total inicial de substrato 
(S, + X,/ Y). 

máxima tasa específica de captación de 
nutrientes: uN (máxima tasa específica 
de captación de nitrógeno) y uP (máxima 
tasa específica de captación de fósforo). 
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1/Y: composición de la biomasa en un instante 
dado. 


Sp: contenido en nitrógeno o fósforo particu- 
lado en un instante f. 
Sp: Concentración inicial de substrato parti- 


culado. 
AY: procesos de acumulación o pérdida de 
reservas en la biomasa. 


S e concentración de sustrato obtenida de 
sustituir la ux (cuadro 1) en la ecuación 
(6). 

S.: concentración de sustrato obtenida de 


us 
sustituir la uS (cuadro 2) en la ecuación 


(6). 
Resultados y discusión 
Luz y temperatura 


La radiación solar y la temperatura se monito- 
rearon a lo largo de todo el experimento, obte- 
niéndose una radiación solar media de 562.87 
+ 145.94 umol m2? d* (p = 95%, n = 107), y una 
temperatura media en el interior y exterior del 
HRAP de 20 + 0.71 *C y de 21 + 0.99 *C (p = 95%, 
n = 107), respectivamente. 


pH 


El pH al inicio del experimento fue de 8.2 y fue 
aumentando hasta que aproximadamente a 
las 72 h de experimento se estabilizó en torno a 
un pH medio de 10.27 + 0.08 (p = 95%, n = 24), 
pH superior al obtenido por Arbib et al. 
(2012) para este mismo fotobiorreactor HRAP 
(pH = 9.32 + 0.30), pero utilizando agua residual 
en lugar de medio sintético, o al valor de 8.5 
obtenido por Ketheesan y Nirmalakhandan 
(2012), cultivando Scenedesmus sp., en medio 
basal Bold 's en un raceway de 23 1 con ilumina- 
ción artificial en laboratorio. Estas diferencias 
pueden deberse a la hora a la que se realizó la 
medición del pH. 

El valor de pH registrado es muy depen- 
diente de la hora a la que se hace la medición. 
En nuestro caso, el pH se medía todos los días 
aproximadamente a la misma hora (9:00-11:00), 
pero el pH evoluciona a lo largo del día, como 
han indicado algunos autores (Ketheesan éz Nir- 
malakhandan, 2012; Putt, Singh, Chinnasamy, éz 
Das, 2011). Para comprobar dicha evolución, se 
realizó un registro completo del pH durante 24 
horas (figura 2). 
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Figura 2. Evolución del pH en 24 horas desde el 30/05 a las 7:00 horas al 31/05 a las 7:00 horas. 
El 30/05, la salida del Sol se produjo a las 7:09 horas y la puesta a las 21:37 horas. 
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Figura 3. Evolución del bloom de microalgas en medio combo. Los rombos hacen referencias a los datos experimentales y la 


línea a los datos predichos por el modelo de Verhulst. 


En la figura 2 se observa perfectamente como 
el pH aumenta durante el día, alcanza su pH 
máximo en la puesta de sol y disminuye de 
modo gradual durante la noche. Estas variacio- 
nes de pH se deben a los procesos biológicos 
de fotosíntesis (consumo de CO,) y respiración 
(generación de CO,) llevados a cabo por las 
microalgas (Margalef, 1983). 


Crecimiento de biomasa 


En la figura 3 se muestra el crecimiento del bloom 
de microalgas a lo largo de todo el experimento. 
Las microalgas crecen conforme a las fases de 
crecimiento microbiano en cultivos en discon- 
tinuo: fase de latencia, fase exponencial, fase 
estacionaria y fase de declive (Prescott, Harley, 
éz Klein, 2004). En la figura 3 se observan con 
claridad las fases de crecimiento exponencial y 
estacionaria, mientras que las fases de latencia 
y declive no son apreciables. La fase de latencia 
probablemente no se aprecie debido a la simili- 
tud con las condiciones ambientales bajo las que 
se obtuvo el inóculo y la fase de declive debido a 
que el experimento se dio por finalizado cuando 
se alcanzó la fase estacionaria, aproximadamen- 
te a las 300 h de haberse iniciado el experimento. 


La cinética de crecimiento de la biomasa se 
ha analizado de acuerdo con el modelo de Ver- 
hulst (Verhulst, 1838) y sus parámetros cinéticos 
aparecen recogidos en el cuadro 1. De acuerdo 
con este modelo, el crecimiento microbiano pue- 
de ser expresado como una curva sinusoidal, 
como describe la ecuación (1), donde X sería la 
concentración de microalgas en un instante t; 
X,, la concentración inicial de microalgas; X,, 
la concentración máxima de microalgas que el 
sistema alcanzaría en batch; u, la máxima tasa de 


crecimiento específico, y t es el tiempo: 


ut 
XX me 


“E “1 


Por otra parte, con el fin de comparar los 
resultados obtenidos con otros experimentos, 
también se ha calculado la productividad 
volumétrica (cuadro 1), que se define como la 
biomasa producida por volumen de reactor y 
unidad de tiempo, es decir, nos muestra la capa- 
cidad que tiene un reactor de producir biomasa 
bajo condiciones de operación específicas. Este 
parámetro puede ser calculado a partir de la 
ecuación (2) (Ruiz et al., 2012a) utilizando los 
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parámetros cinéticos del modelo de Verhulst 


(cuadro 1): 
-(0.9X -1.1X 
Pon o X, 1 E dl 
ln ( most ,) 
11X, 


El modelo cinético de Verhulst se ajusta bien 
a los datos experimentales (figura 3), siendo 
el coeficiente de determinación (R?) de 0.946 
(cuadro 1). Además, el valor de X, medido ex- 
perimentalmente en el laboratorio coincide con 
el predicho por el modelo de Verhulst, siendo 
ambos de 0.072 y 0.073 g SS 1”, respectivamen- 
te, lo que también corrobora el buen ajuste al 
modelo cinético de Verhulst. La Xm es de 0.284 
g SS L”, prácticamente igual a la obtenida por 
Arbib et al. (2013), que fue de unos 0.300 g SS 
1, La similitud entre ambos puede deberse al 
hecho de que en los dos experimentos se uti- 
lizó el mismo fotobiorreactor, aunque Arbib et 
al. (2013) emplearon agua residual pretratada 
y un monocultivo de Scenedesmus obliquus. La 
máxima tasa de crecimiento específico obtenida 
ha sido de 0.442 d”, unas cuatro veces superior 
a la obtenida por Radmann, Reinehr y Costa 
(2007), quienes obtuvieron una u de 0.038 - 0.138 
d* para Spirulina platensis cultivada en medio 
Zarrouk, en un raceway con un volumen de seis 
litros a una concentración de 0.15 g 1* y en con- 
diciones controladas de luz y temperatura. En 
cuanto a la productividad volumétrica (0.025 g 
SS 1" d1), ésta sí es comparable a los 0.028-0.046 
g 1? d obtenidos por Radmann et al. (2007), 
pero inferior a la obtenida por Ketheesan y 
Nirmalakhandan (2012), que fue de 0.085 g 1 d* 
para Scenedesmus sp. cultivada en medio basal 
Bold's en un raceway de 23 litros con ilumina- 
ción artificial en el laboratorio, probablemente 
porque estos autores airearon el reactor con aire 
enriquecido con CO, (1%). 


Por otro lado, a partir de los parámetros ciné- 
ticos de crecimiento obtenidos en el experimento 
en discontinuo, se ha predicho el valor máximo 
de concentración de biomasa (ecuación (3)), y 
las productividades volumétrica (ecuación (4)) 
y superficial (ecuación (5)) del fotobiorreactor 
operando en continuo (Ruiz et al., 2012b): 


x-x, (1-5) (3) 
u0 
a 1, 
px, (0 a) (4) 
U 
P, =P, Buap (5) 


Donde 0 es el tiempo de retención hidráulico 
(TRH) al que opera el fotobiorreactor y h,yz,p 
es la profundidad de la lámina de agua del 
fotobiorreactor. Considerando que la máxima 
productividad en continuo se alcanza para un 
TRH de 2/u (Ruiz et al., 2012b). 

El valor de concentración de biomasa (ecua- 
ción (3)) y máxima productividad volumétrica 
(ecuación (4)), operando en continuo a un TRH 
de 4.54 d (2/1) serían de 0.142 y 0.031 g SS 1" 
d*, respectivamente. El TRH óptimo al que 
se tiene la máxima producción de biomasa es 
aproximadamente la mitad al obtenido por Ar- 
bib et al. (2012) para este mismo fotobiorreactor 
HRAP, pero utilizando agua residual en vez 
de medio sintético. Son muchos los autores 
que consideran que el TRH óptimo al que se 
produce la máxima producción de biomasa se 
encuentra entre los 8 y 10 d (De-Godos, Blanco, 
García-Encina, Becares, £ Muñoz, 2009; García 
et al., 2006). El hecho de que nuestro TRH óp- 
timo al que se produce la máxima producción 
de biomasa sea inferior se debe al elevado valor 
de u (cuadro 1), como consecuencia del uso de 
medio sintético como forma de cultivo. 


Cuadro 1. Parámetros cinéticos para el crecimiento de un bloom de microalgas según el modelo de Verhulst. 


X, (8 SS 1”) X, (858 1") 


u (d”) 


R? 


P, (8 SS 1" d) 


0.073 0.284 


0.442 


0.025 0.946 
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El caudal al que debería ser alimentado el 
fotobiorreactor para obtener la máxima pro- 
ducción de biomasa funcionando en continuo 
es de 99 1 d* y la concentración de biomasa en 
el efluente sería de 141.75 mg 1”. Esta concen- 
tración de biomasa en el efluente es también 
aproximadamente la mitad a los 250-300 mg 
11 que Tredeci (2004) determinó como concen- 
tración de biomasa típica para el efluente de 
un fotobiorreactor HRAP. Pero se debe tener 
en cuenta que esta producción de biomasa es 
para raceway a escala industrial, en los que el 
efecto sombreado creado por las paredes del 
fotobiorreactor es mínimo y, en nuestro caso, se 
trata de un raceway piloto donde el efecto de 
sombreado en los canales debido a su estrechez 
tiene como consecuencia una disminución de 
productividad, lo que pone en evidencia el 
problema de escalado que presentan los reacto- 
res tipo HRAP, al igual que han indicado otros 
autores (Ketheesan €: Nirmalakhandan, 2011). 

La productividad superficial estimada es de 
9.4 g m? d” (ecuación (5)), aproximadamente un 
30% más baja de la productividad superficial 
típica para un HRAB que es de 12-13 g m? d* 
(Tredeci, 2004). De-Godos et al. (2009) operaron 
dos HRAP para tratar purines diluidos a un 
TRH de 10 d, el doble al estimado para este 
experimento, y obtuvieron una productividad 
superficial de 21 a 28 g m? d”. García et al. 
(2006) obtuvieron una productividad superficial 
media mensual de 12.7 g m” d* para un HRAP 
operando a un TRH de 10 d y 14.8 g m? d* 
para otro HRAP, pero operando a un TRH de 
8 d, mientras que Arbib et al. (2012) obtuvieron 
una productividad superficial de 8.26 + 1.43 g 
m? d*, operando a un TRH de 8 d el mismo 
fotobiorreactor HRAP que el utilizado para este 
experimento, resultando que sí se encuentra en 
la misma línea que los estimados en el presente 
estudio, lo cual corrobora que la baja produc- 
tividad se debe a aspectos relacionados con la 
escala del fotobiorreactor. 

En este estudio también se estimó la eficacia 
fotosintética del sistema operando en disconti- 
nuo, que fue de 0.650%, y la que habría tenido 
si hubiese operado en continuo a un TRH de 4.5 
d, esto es, sólo un 0.845%. Ésta se estimó a partir 


de la radiación PAR media medida, partiendo 
de que la energía solar en el intervalo de 400 a 
700 nm es de un 48% de la energía solar total 
(Zhu, Long y Ort, 2008) y que el poder calorífico 
superior medio de las microalgas es de 20 MJ K* 
(Illman, Scragg, € Shales, 2000). La eficiencia 
fotosintética típica de un HRAP es de 1.3 a 2.4% 
(Park, Craggs, €: Shilton, 2010), superior a la 
obtenida en este estudio, lo que puede deberse, 
de nuevo, al efecto sombreado creado por las 
paredes del fotobiorreactor. 


Eliminación de nutrientes 


El consumo de nitrógeno y fósforo por las 
microalgas se ha analizado de acuerdo con el 
modelo cinético de Romero de consumo de 
substratos por microorganismos en reactores en 
discontinuo (Quiroga et al., 1999). De acuerdo 
con este modelo, Ruiz et al. (2012a) desarrolla- 
ron un modelo para predecir la captación de 
nutrientes por microalgas en un fotobiorreactor 
operando en discontinuo, modelo PhBT, en el 
que la captación de nutrientes puede expresarse 
a partir de la ecuación (6), donde S, es la con- 
centración de substrato inicial (nutrientes en 
el medio de cultivo); S, es la concentración de 
substrato no asimilable; St, es la concentración 
total inicial de substrato (nutrientes en el medio 
de cultivo (S,) y nutrientes en el inóculo (X,/ Y) 
(ecuación (7)); 1 es la máxima tasa específica de 
captación de nutrientes, y t es el tiempo: 


S (E +S Je 
10) Y 10) 
- (6) 


(S, -S,,)- s - E +5, Je 


S o (7) 
Y 


Para el cálculo de los parámetros cinéticos 
S,, S,, y u se ha partido del valor conocido de 
S,, que se puede estimar a partir de la ecuación 
(7) con los datos experimentales obtenidos del 


na 


análisis elemental de la biomasa empleada como 
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inóculo (1 / Y, = 0.062 mg N mg SS” y 0.009 mg P La capacidad de eliminación de nitrógeno y 
mg SS”) y el valor de la concentración inicial de fósforo por el bloom de microalgas se muestran 
biomasa obtenido de la modelización cinética en las figuras 4 y 5, y sus parámetros cinéticos 
del crecimiento, X, =73 mg 1* (cuadro 1). aparecen recogidos en el cuadro 2. 


0? za E ES ; an F = > 
0 50 100 150 200 250 300 350 400 


Tiempo (h) 


Figura 4. Eliminación del nitrógeno total soluble del medio en el HRAP por el bloom de microalgas cultivadas en combo. 
Los rombos hacen referencias a los datos experimentales y la línea a los datos predichos por el modelo cinético. 


2.57 


0 50 100 150 200 250 300 
Tiempo (h) 


Figura 5. Eliminación del fósforo total soluble del medio en el HRAP por el bloom de microalgas cultivadas en combo. 
Los rombos hacen referencias a los datos experimentales y la línea a los datos predichos por el modelo cinético. 


Cuadro 2. Parámetros cinéticos de eliminación del nitrógeno y fósforo total mediante el modelo PhBT. 


S, (mg 1) s, (mgl) S, (mg 1) u (di) R? 
N 14.101 2.615 18.624 0.359 0.934 


ll 1.914 0.046 2.666 1.313 0.982 
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Eliminación del nitrógeno total 


El modelo cinético se ajusta bien a los datos ex- 
perimentales (figura 4), siendo el coeficiente de 
determinación R? de 0.934 (cuadro 2). Además, 
la S, medida de forma experimental en el labo- 
ratorio se aproxima bastante con la predicha 
por el modelo, siendo ambas de 16.392 y 14.101 
mg N,, L”*. La concentración de S, , que es el 
contenido de nitrógeno total soluble que tiene 
el efluente a la salida del fotobiorreactor, es de 
2.615 mg N, 1”, habiéndose reducido la con- 
centración de nitrógeno total en un 81.5% a la 
salida del efluente y siendo aproximadamente 
un 75% inferior al límite más restrictivo de ver- 
tido de 10 mg N,, 1* que establece la Directiva 
98 /15/CE para aguas residuales urbanas. La u,, 
es de 0.359 d”, inferior a los 1.4-2.5 d obtenidos 
por Ruiz et al. (2012a), pero hay que tener en 
cuenta que Ruiz et al. (2012a) cultivaron Chlore- 
lla vulgaris en diferentes medios de cultivos en 
frascos Pyrex de borosilicato de dos litros, y en 
condiciones controladas de luz y temperatura 
en el laboratorio. 


Eliminación del fósforo total 


El modelo cinético se ajusta bien a los datos 
experimentales (figura 5) y coeficiente de 
determinación R? de 0.982 (cuadro 2); ade- 
más, la S, medida de forma experimental en 
el laboratorio es muy parecida a la predicha 
por el modelo, siendo ambas de 2.122 y 1.914 
mg P, 1”, respectivamente. La concentración 
de S, es de 0.046 mg P, 1*, reduciéndose la 
concentración de fósforo prácticamente en 
su totalidad a la salida del efluente y siendo 
inferior al valor más restrictivo de concentra- 
ción en vertido de 1 mg P,, 1* que establece 
la Directiva 98/15/CE. La u, es de 1.313 di, 
mayor a la u, (cuadro 2) obtenida para el 
nitrógeno, lo cual significa que la velocidad 
específica de consumo de fósforo por las 
microalgas es mayor. La u,, al igual que en el 
caso del nitrógeno, es muy inferior a los 2-8.7 
d obtenidos por Ruiz et al. (2012a). 


Composición de la biomasa 
Evolución del nitrógeno en la biomasa 


Según el modelo PhBT (Ruiz et al., 2012a), el 
valor de u del modelo de crecimiento de bio- 
masa y de consumo de nutrientes son el mismo 
parámetro; pero en nuestro caso, al compararlos 
se observa cómo no coinciden, ya que la u del 
modelo de crecimiento de biomasa es de 0.441 
d* (cuadro 1) y la de consumo de nitrógeno es 
de 0.359 d* (cuadro 2). 

En la figura 6 se ha representado la evolución 
de los nutrientes predicha por la ecuación (6) y 
usando los valores de S, y $, recogidos en el 
cuadro 2, pero empleando en un caso el valor 
de u obtenido de la modelización de la biomasa 
(u,) y en otro el valor fruto de la modelización 
del consumo de nitrógeno (u,). 

En la figura 6 se observa cómo parece ser que 
las microalgas consumen nitrógeno (línea negra) 
a una velocidad superior a la que realmente está 
descendiendo la concentración de nitrógeno en 
el medio (línea gris). Pero esto no es posible, 
pudiéndose deber a que la concentración de 
nitrógeno en la biomasa (1/Y,) varía con el 
tiempo (ecuación (8)): 


dX _y 4N 6 
dt dt 

En la figura 8 se puede observar cómo, en 
efecto, los valores experimentales de contenido 
en nitrógeno de la biomasa no son constantes 
en el tiempo. 

En el presente trabajo se ha desarrollado una 
ecuación que justifique los procesos de acumu- 
lación o pérdida de nutrientes en la biomasa a 
lo largo del proceso de fotobiodepuración. Para 
ello se parte de la ecuación (9), en la que se ex- 
presa el valor de la composición de la biomasa 
en un instante dado (1/Y) como: 

e (9) 
Y Xx 

Donde X es la concentración de la biomasa 
en un instante t y S, es el contenido en nitrógeno 
o fósforo particulado en un instante f. El conte- 
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0+ 7 7 7 Y 7 7 1 
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Tiempo (h) 
a.  Npred (uN) eu N pred (uXx) 


Figura 6. Predicción de la evolución del nitrógeno según el PhBT con la máxima tasa de crecimiento específico de la biomasa 
(u,) y con la máxima tasa específica de captación de nitrógeno (u,). 


nido en nitrógeno o fósforo particulado puede 


expresarse como S,= S,, + (S, - 5), donde S,, es ; E Jus -S,s) 
la concentración inicial de substrato particulado; =SA—9 
S, es la concentración inicial de substrato solu- Y; Xx (12) 
ble, y S es la concentración de substrato soluble 
en un instante t. Por lo tanto, sustituyendo en 
la ecuación (9) obtenemos: : E Jus -S.x) 
13 
(%)+(5,-5 Y A e 
1 1% 
Y X (10) Reordenando términos, obtenemos que: 
El problema es que el valor de S no es el 
mismo si empleamos la velocidad específica a+) E 2.(5 $=8 2) (14) 
de crecimiento derivada de la modelización Y) X*" j ” 
de la biomasa (u,) que si empleamos el valor 5 
derivado del consumo de substrato (u, O 41), y Sustituyendo la ecuación (6) en la ecuación E 
es precisamente esa diferencia la que define los (14) y reordenando términos, tenemos la ecua- El 
procesos de acumulación o pérdida de reservas ción (15), donde A=S,-S, y B=-X,/Y; A 
en la biomasa durante el proceso (AY) (ecuación 0 
(11) E 
197 ust ] 
a an a lo a ) | : 
YES Ye Ye A -A-Be[e + el )+B ee gus! y 
5 
Donde: (15) 3 
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En la figura 7 puede verse cómo los datos que 
predice el modelo (ecuación (10)) no coinciden 
con los datos experimentales de composición de 
la biomasa, siendo el contenido experimental de 
nitrógeno en la biomasa menor al esperado, lo 
que puede deberse a que no todo el nitrógeno 
eliminado en el agua residual es consumido por 
las microalgas, sino que existen otros procesos 


que también participan en la eliminación del 
nitrógeno. 

Para cuantificar la cantidad de nitrógeno que 
no ha sido consumido por las microalgas y que 
ha sido eliminada del medio, se determinará el 
valor de y, necesaria para que se dé un buen 
ajuste de los datos experimentales de composi- 
ción (cuadro 3; figura 8). 


1/Yx 0) 


Ed sd 1 / Y experimental 


su. 1 / Y predicho 


50 100 150 200 
Tiempo (h) 


Figura 7. Contenido de nitrógeno en la biomasa con respecto al tiempo, comparando los valores obtenidos de manera 


experimental con los predichos por este estudio. 


4 1/Y experimental 


===] / Y predicho 


0 50 100 150 200 
Tiempo (h) 


250 300 350 


Figura 8. Contenido de nitrógeno en la biomasa correctamente ajustado con base en la evolución de la composición 


de la biomasa. 


Cuadro 3. Parámetros cinéticos de eliminación del nitrógeno ajustados a los datos experimentales de composición de biomasa. 


gua, vol. VI, núm. 3, mayo-junio de 2015 


S, (mg N, 13) S, (mg N, 1) S, (mg N, 13) 


uy (d=) 1/Y (mg SS mg N”) R? 


18.494 


13.968 9.032 


0.356 0.062 0.995 
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Según estos, el consumo del nitrógeno 
soluble en el medio debería haber sido el 
representado en la figura 9 (línea negra), por 
lo que parece ser que existe otro proceso com- 
plementario al consumo biológico de nitrógeno 
por las microalgas que también participa en 
la eliminación del nitrógeno del medio. Estos 
procesos complementarios podrían ser bióticos, 
como la desnitrificación, o abióticos, como el 
stripping. El stripping es un proceso abiótico en 
el que a temperatura ambiente y pH constante 
superiores a 8, como sería nuestro caso, el NH,' 
se transforma en amoniaco gaseoso (NH,) y 
se libera a la atmósfera (Margalef, 1983). Pero 
esta hipótesis queda descartada, pues la forma 
mayoritaria de nitrógeno en el medio era NO),, 
siendo la concentración inicial de NO, y NH,* 
en el medio de 17.4 y 0.625 mg l”, respectiva- 
mente. La desnitrificación es la transformación 
biológica del nitrato en gas nitrógeno (N,), óxido 
nítrico (NO) y óxido nitroso (N,O). Todos estos 
compuestos son gaseosos y se liberan normal- 
mente a la atmósfera. Pero éstas son transforma- 
ciones que se producen en condiciones anóxicas 
(Margalef, 1983) y, en nuestro caso, el sistema es 
oxidante en extremo, dada la elevada actividad 
fotosintética en el mismo. Descartadas estas dos 


opciones, queda otro proceso biológico, menos 
conocido, aunque no por ello menos importante, 
como es el hecho de que los cultivos de algas 
son fuentes de emisión de CH, y N,O. De ma- 
nera reciente, autores como Harter et al. (2013), 
o Ferrón, Ho, Johnson y Huntley (2012) han 
cuantificado los flujos de N,O desde el agua al 
aire. Ferrón et al. (2012), cultivando Staurosira sp. 
en un raceway de 64 000 litros expuestos a luz 
solar y a la atmósfera, determinaron que para 
raceway fertilizados durante todo el periodo de 
cultivo, como es nuestro caso, el flujo de N,O del 
agua a la atmósfera sería de 342.86 umoles m? 
d”, y Harter et al. (2013), cultivando Dunaliella 
salina en un reactor de 42 litros, suministrándole 
gas de combustión en condiciones controladas 
de luz y temperatura, determinaron que este 
flujo era de 2 711 ug NI” d*”. 

En la figura 10 se presenta la evolución del 
nitrógeno en el medio debido a los procesos 
complementarios, obtenido por diferencia de 
las dos curvas anteriores (figura 9). Se puede 
concluir que del total del nitrógeno eliminado, 
un 43.5% ha sido consumido por las microalgas, 
mientras que el 56.5% restante fue eliminado 
por procesos complementarios; posiblemente 
parte de este nitrógeno se emitió a la atmósfera 


su. N biótico 


— == N total 


50 100 150 


250 


Tiempo (h) 


Figura 9. Eliminación del nitrógeno en el medio, comparando los valores predichos con base en los datos experimentales 
de evolución del nitrógeno soluble en el medio de cultivo (en gris) (N total), y los predichos con base en la evolución de la 
composición de la biomasa (en negro) (N biótico). 
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Figura 10. Evolución del nitrógeno en el medio debido a los procesos complementarios. 


en forma de N,O. Si suponemos que todo el 
nitrógeno que no ha sido consumido por las 
microalgas se ha emitido a la atmósfera en 
forma de NO, se obtiene que este flujo ha sido 
de 5 200.16 umol N,O m” d*, que es 15 veces 
superior a lo obtenido por Ferrón et al. (2012), 
o lo que es lo mismo, 485.35 1g N 1* d*, que es 
cinco veces inferior a lo que obtuvo Harter et 
al. (2013). 


Evolución del fósforo en la biomasa 


En el caso del fósforo, la u del modelo de creci- 
miento de biomasa y la de consumo de fósforo 
tampoco coinciden, siendo ambas de 0.44 d* 
(cuadro 1) y de 1.31 d* (cuadro 2), respectiva- 
mente. 

En la figura 11 se representa la evolución 
del fósforo predicha a partir de la u del modelo 
de consumo de nutrientes y del crecimiento 
de biomasa (ecuación (6)), y a diferencia del 
nitrógeno, puede observarse cómo la biomasa 
consume fósforo (u,) más lentamente de lo que 
lo realmente lo hace el sistema (u,). El hecho de 
que ambas curvas no coincidan es posible que 
se deba al efecto combinado de dos procesos: 
uno biológico, conocido como “luxury uptake”, 
esto es, que las microalgas tienden a consumir 


más nutrientes del que necesitan para crecer 
acumulándolos en forma de reserva (Ruiz el al., 
2012a); pero también puede ser que participen 
procesos abióticos, como la precipitación de 
fosfatos, dado el pH alcanzado en los reactores. 

En la figura 12, al igual que para el nitró- 
geno, también puede verse cómo los valores 
experimentales de contenido en fósforo de la 
biomasa varían con el tiempo (ecuación (8). 
En esta figura también se han representado los 
valores de 1/Y, predichos por la pérdida de 
fósforo en la biomasa a partir de la ecuación 
(10), y se observa cómo los datos que predice el 
modelo se ajustan a los datos experimentales de 
composición de la biomasa, con excepción del 
valor máximo observado de alrededor de 1.9% 
al cabo de 22 h. Este resultado podría indicar 
que la rápida eliminación observada para el 
fósforo en las primeras 70 h (figura 5) podría 
deberse a un proceso de precipitación, lo cual 
haría que la composición en fósforo del material 
particulado del reactor (1.9%) sea muy superior 
al valor máximo típico en microalgas, que no 
suele superar el 1%. De manera posterior, este 
fósforo extracelular se iría asimilando por parte 
de las microalgas, para finalmente formar parte 
del material celular en una importante propor- 
ción en los estadios finales del crecimiento en 
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Figura 11. Predicción de la evolución del fósforo según el PhBT con la máxima tasa de crecimiento específico de la biomasa (u,) 


y con la máxima tasa específica de captación de fósforo (up). 
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Figura 12. Contenido de fósforo en la biomasa con respecto al tiempo, comparando los valores obtenidos de manera 
experimental con los predichos por este estudio. 


batch, coincidiendo los valores experimentales 
con los predichos por el modelo. Esta hipótesis 
debería ser confirmada mediante la compa- 
ración de estos resultados con los obtenidos 
del análisis elemental de fósforo a la biomasa, 
previo lavado de la misma con una disolución 
ligeramente ácida, para discernir entre las dos 
formas de fósforo en la biomasa. 


Conclusiones 


El cultivo de un bloom de microalgas en un fo- 
tobiorreactor piloto tipo HRAP permite reducir 
en más de un 80% las concentraciones de nitró- 
geno y fósforo presentes en el medio sintético, 
lo que apunta a que se trata de una tecnología 
con elevada potencialidad como tratamiento 
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terciario para la eliminación de nutrientes en 
estaciones depuradoras de aguas residuales. 
El modelo cinético PhBT es una herramienta útil 
para predecir la cinética de crecimiento de la 
biomasa, y de consumo de nitrógeno y fósforo 
en fotobiorreactores operando en discontinuo. 
Se ha podido comprobar de forma experimental 
que el modelo PhBT además predice con bastan- 
te aproximación la evolución de la composición 
de la biomasa en términos de contenido en 
nitrógeno y fósforo, lo cual puede resultar muy 
útil a la hora de identificar otros procesos que 
compiten con los biológicos en la eliminación de 
nutrientes en los fotobiorreactores. Por último, 
el conocimiento de la evolución de la riqueza de 
nitrógeno y fósforo en la biomasa es un aspecto 
clave a la hora de establecer estrategias de ope- 
ración orientadas hacia la minimización en la 
producción de biomasa algal, manteniendo una 
elevada velocidad de consumo de nutrientes. 
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